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1. Introduction 

We present the recent accurate determination of the average up/down, strange and charm 
quark masses performed by ETMC with Nf = 2 maximally twisted mass Wilson fermions. The 
high precision of this analysis is mainly due to the extrapolation of the lattice results to the contin- 
uum limit, based on data at four values of the lattice spacing (a ~ 0.098, 0.085, 0.067, 0.054 fm), to 
the well controlled chiral extrapolation, which uses simulated pion masses down to M % ~ 270 MeV, 
and to the use of the non-perturbative renormalization constants calculated in [^]. The only system- 
atic uncertainty which is not accounted for is the one due to the missing strange and charm quark 
vacuum polarization effects. However, a comparison of Nf = 2 results for the up/down and strange 
quark masses to already existing results from Nf = 2 + 1 quark flavor simulations [Q] indicates 
that, for these observables, the error due to the partial quenching of the strange quark is smaller 
at present than other systematic uncertainties. The same conclusion is expected to be valid for the 
effects of the strange and charm partial quenching in the determination of the charm quark mass. 
In this respect we mention that simulations with Nf = 2 + 1 + 1 dynamical flavors are already being 
performed by ETMC and preliminary results for several flavor physics observables have been re- 
cently presented []5j, ^. For more details on the ensembles of Nf = 2 gauge configurations used in 
the analysis and the values of the simulated light, strange and charm quark masses we refer to 

The calculation of the averaged up/down quark mass, based on the study of the pion mass and 
decay constant, has closely followed the strategy of At variance with the latter, however, here 
and in [Q] data at four values of the lattice spacing have been used. For the strange quark mass, the 
main improvement with respect to our previous work which used data at a single lattice spacing 
only, is the continuum limit. Moreover, the chiral extrapolation has been performed by using either 
SU(2)- or SU(3)-Chiral Perturbation Theory (ChPT). In order to extract the strange quark mass 
we have used both the kaon mass and the mass of the (unphysical) T] s meson composed of two 
degenerate valence strange quarks. For the charm quark mass, similarly to the strange quark, we 
have used several experimental inputs to extract its value: the mass of the D, D s and T] c mesons. 

The results that we have obtained for the quark masses are, in the MS scheme, 

7n Hd (2GeV) = 3.6(2) MeV, 
m s {2 GeV) = 95(6) MeV, 

m c (m c ) = 1.28(4) GeV. (1.1) 

We have also obtained for the ratios of quark masses the values 

m s /niud = 27.3(9) and m c /m s = 12.0(3) , (1.2) 

which are independent of both the renormalization scheme and scale. 

Finally, we take the opportunity of these proceedings to present the updated result for the 
bottom quark mass which has been obtained employing the method discussed in [|l|]. The updated 
value is 

mb(Mb) =4.3(2) GeV . (1.3) 
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2. Up/down quark mass 

We have studied the dependence of the pion mass and decay constant on the renormalized 
quark mass, by using the predictions based on NLO ChPT and the Symanzik expansion up to 
&(a 2 ). 1 They include discretization terms of tff(a 2 mi\og(mi)), which receive a contribution from 
the 0(a 2 ) splitting between the neutral and charged pion mass []|] occurring with twisted mass 
fermions. The impact of this correction on the final result for the light quark mass is at the level of 
the fitting error. 

Lattice results for pion masses and decay constants have been corrected for finite size effects 
(FSE) evaluated using the resummed Liischer formulae. The effect of the &(a 2 ) isospin breaking 



has been taken into account also in these corrections QlOp . On our pion data, FSE vary between 
0.2% and 2%, depending on the simulated mass and volume. The inclusion of the pion mass 
splitting in the FSE induces an effect at the level of one third of the statistical error for our lightest 
pion mass at j8 = 3.9 on the smaller volume, and even smaller in the other cases. 

The value of the physical up/down quark mass is extracted from the ratio m^/f 2 using as an 
input the experimental value of the latter ratio. In order to estimate the systematic uncertainty due 
to discretization effects we have performed both a fit without the logarithmic discretization terms, 
and a fit without all &{a 2 ) corrections. Both these ansatze turn out be compatible with the lattice 
data. We find that the result for the up/down quark mass decreases by approximately 2% and in- 
creases of about 6% in the two cases respectively, so that we estimate an overall uncertainty due to 
residual discretization effects of ±4%. For estimating the systematic uncertainty due to the chiral 
extrapolation we have also considered a fit including a NNLO local contribution proportional to 
the light quark mass square. In this case we are not able to determine all the fitting parameters 
and we are thus forced to introduce, on the additional LECs, priors as in 01. In this way we find 
that the result for m u d increases by 6%. We have also included in the final result a 2% systematic 
uncertainty coming from the perturbative conversion of the quark mass renormalization constant 
from the RI-MOM to the MS scheme. This uncertainty has been conservatively estimated by as- 
suming the unknown ^(ot s 4 ) term to be as large as the ^(ot 4 3 ) one, evaluated at the renormalization 
scale ~ 3 GeV, which is the typical scale of the non-perturbative RI-MOM calculation in our 
simulation [S]. Adding in quadrature the three systematic errors discussed above we have obtained 

m ud {2 GeV) = 3.6(1)(2) MeV = 3.6(2) MeV. (2.1) 



3. Strange quark mass 

In this section, we first present the determination of the strange quark mass based on the study 
of the kaon meson mass and then the alternative method based on the study of the r\ s meson. 

In order to better discriminate the strange quark mass dependence of the kaon masses on other 
dependencies, we have firstly interpolated, by using quadratic splines, the lattice data to three 
reference values of the strange quark mass, chosen to be equal at the four lattice spacings: = 
{80, 95 , 110} MeV. Then, at fixed reference strange mass, we have simultaneously studied the 
kaon mass dependence on the up/down quark mass and on discretization effects, thus performing 

'Here and in the following sections we refer to [0] for the expressions of the fitting functions. 
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Figure 1: Left: Dependence of nv\ and m^ s on the renormalized light quark mass, for a fixed reference 
strange quark mass (m™ = 95 MeV ) and at the four lattice spacings. Right: Dependence of m\ and m^ s on 
the squared lattice spacing, form^ = 95 MeV and at the physical up/down mass. 
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Figure 2: Dependence ofm\ and m^, in the continuum limit and at the physical up/down mass, on the 
strange quark mass. The strange mass results are also shown (empty diamonds). 



a combined chiral and continuum extrapolation. In this step, we have considered chiral fits based 
either on SU(2)-ChPT [JTTJ, [f|] or partially quenched SU(3)-ChPT Q. Finally, we have studied the 
kaon mass dependence on the strange quark mass, and determined the value of the physical strange 
quark mass using the experimental value of mk ■ In fig. [I] we show the combined chiral/continuum 
fit based on SU(2)-ChPT, for a fixed reference value of the strange quark mass, as a function of 
the light quark mass (left) and of the squared lattice spacing (right). In fig. || the dependence on 
the strange quark mass is shown, for the SU(2) analysis. The dependencies are shown for the kaon 
squared mass as well as for the rj s squared mass discussed hereafter. 

As an alternative way to determine the strange quark mass we have studied the dependence 



on m s of a meson made up of two strange valence quarks [14]. The advantage of this approach is 
that the mass of this unphysical meson, denoted as T] s , is only sensitive to the up/down quark mass 
through sea quark effects, and thus requires only a very smooth chiral extrapolation. The price to 
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pay is the need for an additional chiral fit to determine the rj s mass at the physical point. In order 
to relate the mass of the t] s meson to the physically observable m n and tuk , we have studied its 
dependence on the kaon and pion masses for different values of the simulated light and strange 
quark masses. We investigated functional forms based on either SU(2)- or SU(3)-ChPT. The two 
fits yield very close results for the r] s meson mass and we quote as final estimate m rjs = 690(3) 
MeV (to be compared to the LO SU(3) prediction m ns = (2m| - m 2 n ) l l 2 = 686 MeV). 

Once the mass of the r\ s meson has been determined, the strange quark mass can be extracted 
by following the very same procedure described for the case of the kaon mass. 

The difference between the determinations based on the K and tj s mesons is about 3%. The re- 
sults obtained from either the SU(2) or the SU(3) fits are practically the same in the analysis based 
on the r\ s and differ by approximately 3% in the kaon case. In order to evaluate the uncertainty 
of the continuum extrapolation we have excluded from these fits the data from the coarser lattice, 
finding a variation of the results of approximately 2%, with the fitting error approximately un- 
changed. The different fits considered for the determination of the up/down mass and of the lattice 
spacing affect the determination of the strange mass at the level of 3%. Finally, we have included 
an uncertainty of 2% related to the truncation of the perturbative expansion in the conversion from 
the RI-MOM to the MS scheme. Combining all these uncertainties in quadrature, we quote as our 
final estimate of the strange quark mass in the MS scheme 

Jn s (2 GeV) =95(2)(6) MeV = 95(6) MeV. (3.1) 

Using the determinations of both the strange and light quark masses, we have also obtained a 
prediction for the ratio m s /m M i, which is both a scheme and scale independent quantity: 

m s /m ud = 27.3(5)(7) = 27.3(9) . (3.2) 

4. Charm quark mass 

The determination of the charm quark mass follows, quite closely, the strategy adopted in the 
determination of the strange quark mass discussed in the previous section. In this case, we have 
used as experimental input the masses of the D, D s and r\ c mesons. 

As for the strange quark case, we have first used a quadratic spline fit to interpolate the 
data at three reference values of the charm mass equal at the four j8 values: JfQ (2 GeV) = 
{1.08, 1.16, 1.24} GeV. In order to fit the meson masses we have considered (phenomenologi- 
cal) polynomial fits, which turn out to describe well the dependence on the light and strange quark 
masses and on the lattice cutoff of the D, D s and r\ c meson masses, at fixed (reference) charm mass 
m c . Then, the value of the physical charm quark mass has been extracted by fitting these data as a 
function of the charm quark mass and using as an input the experimental value of the corresponding 
charmed meson m e ^ p = 1.870 GeV, m^ p = 1.969 GeV, m e £ p = 2.981 GeV. For the charm mass 
dependence, a constant plus either a linear or a 1 jm c term have been considered for describing data 
of the D, D s and T] c mesons. Both choices are found to describe very well the lattice data. In fig. || 
we show the dependence of the D, D s and f] c masses on the light quark mass at a fixed reference 
charm mass, for the four /3's. For the D s and r\ c mesons, which contain the light quark in the sea 
only, this dependence is not significant within the statistical errors. In fig. ^| (left) the meson masses 
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Figure 3: Left: Dependence of mo (left) and mj) s and m-q c (tight) on the light quark mass, at fixed reference 
charm quark mass (m£ =1.16 GeV) and for the four simulated lattice spacings. For the D s meson the 
strange quark mass is fixed to the reference value m'f^ = 95 MeV. 




Figure 4: Left: Dependence of mo, tno s and m T)c , at fixed reference charm quark mass (m c J = 1.16 GeV) 
and at physical up/down and strange quark mass, on the squared lattice spacing. Right: Dependence of mo, 
mo s and m^, in the continuum limit and at physical up/down and strange quarks, on the charm quark mass. 
The charm mass results from the three determinations are also shown (empty diamonds). 

at physical light and strange quark masses are shown as a function of a , for a reference value of 
the charm quark mass. Finally, fig. |] (right) shows the dependence of the D, D s and r\ c masses on 
the charm mass and the interpolation to the physical charm. 

In order to evaluate the systematic uncertainty, we have summed in quadrature the approxi- 
mately 1% spread among the three determinations from the D, D s and T] c mesons, the 1.5% uncer- 
tainty due to discretization effects (estimated by excluding the data from the coarser lattice) and 
the 2% uncertainty coming from the perturbative conversion of the renormalization constants from 
the RI-MOM to the MS scheme. We quote as our final result for the charm quark mass in the MS 
scheme 

m c (2 GeV) = 1.14(3)(3) GeV = 1.14(4) GeV -> m c {m c ) = 1.28(4) GeV, (4.1) 
where the evolution to the more conventional scale given by m c itself has been performed at 
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N 3 LO [15] with Nf = 2, consistently with our non-perturbative evaluation of the renormalization 
constant. Our result is in good agreement with the HPQCD result m c (m c ) = 1.268(9) GeV [juj], 
with a larger uncertainty in our determination, and with the recent sum rules determination m c im c ) = 



1.279(13) GeV of [|17|]. We have also provided a prediction for the scheme and scale independent 
ratio 

m c /m s = 12.0(3) . (4.2) 

5. Bottom quark mass 

In these proceedings we also take the opportunity to provide an updated value for the b-quark 
mass following the method presented in [j|]. This method consists in the calculation of the b-quark 
mass using suitable ratios of the heavy-light pseudoscalar meson masses which, by construction, 
have an exactly known infinite mass limit. In the present study we have employed simulation data at 
three j8 values as in [gl, namely 3.80, 3.90 and 4.05, but using now the final value of the quark mass 
renormalization constants [[J and the same statistics as for the other quark mass determinations. In 
this way we have obtained the updated estimate 

mb{fhb) = 4.3(2) GeV . (5.1) 

This value is about la smaller, with an improved accuracy with respect to the result given in [[I]]. 
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